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We analysed the binomial multiplicity moments of the neutral pi- 
ons, using an extension of the generating functional technique for de- 
tection losses. We applied this model-independent method to the indi- 
vidual 7 weights of 10000 events of tt~p interactions at 40 GeV/c. We 
compared the obtained results to those of 250 GeV/c. We used the 
FRITIOF and a shifted KW distribution to describe the data. 

Introduction 

Investigation of multiplicity distributions have so far been done mostly 
for charged particles. A comprehensive study of functional forms and 
fits for data have been lately reviewed in Warsaw Less informa- 
tion is available on production of 7r° meson. The 7r° decay product, 
gammas, may be observed in bubble chamber with low efficiency. 

We shall analize the moments of multiplicity distributions of 7-s 
for Tr~p and TT~n interactions at 40 GeV. We will use the data from 
the Dubna 2m propane bubble chamber. The statistics includes about 
10000 events for tt~p and 3600 events for TT~n interactions. We have 
25% mean efficiency Q . Every individual conversion weight of 7 is at 
our disposal from the data summary tape (DST). 
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Since the detection probability is lower than 100% the measured 
distribution is different from the true one. The problem is the following: 
how to take into account this difference in the analysis of the data. 



1 General method of the correction of de- 
tector losses. 

The generating functional technique is both extremely general and use- 
ful. On the one hand it can be used to prove important theorems [3|, 
on the other hand it permits the description of detection losses too [ 1 1 . 
In order to give some insight into this problem we will show a general 
model independent method by Diosi We shall recall some state- 
ments from these papers. 

The true n-particle exclusive distribution s^"^ with the proper nor- 
malization is the following: 



(fci, ...,k^)dki,... dkn = n!p„ (1) 



where p„ is the probability of fixed n multiplicity and A:„ is the mo- 
mentum of the n-th particle. 

With the aid of the detection probabilities lo — uj{k)-s we can de- 
scribe the measured exclusive distribution 

s(")(A:i,...,fc„) = V , \, / g<"^(fci,...,fc„)L^(fci)...c^(fc„)- 
■ui{kn+i)dkn+i ■ . . U}{km)dkm (2) 

where uj = 1 — lu. 

By definition the generating functional: 

F [h{.)] ^Y.^J (^1' • ■ ■ ' Hki)dki . . . h{k^)dkn (3) 

The exclusive distribution can be expressed by the derivatives of the 
generating functional: 



^ (k\ , . . . , k„ 



' 5h{ki)...5h{kn) 
The generating functional of the measurable distribution: 



(4) 



F [K-)] ^Y.^J ^^"^ (^1' ■ ■ ■ ' Hki)dki . . . h{kn)dkn (5) 
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Using these eqs. the fundamental reconstruction formula can be ob- 
tained as: 

F[h{.)]^F[cui.)h{.) + l-iu{.)] (6) 

F[hi.)]^F[u;-\.)h{-) + l-^-\-)] (7) 

We can check these formulae: if the argument h{.) or h{.) = 1 then 
= F[l]for arbitrary uj{.), and if uj{.) = 1 then F = F. We should 
note that and (0) is a generalization of Nifenecker's results for the 
generating function if h{.) z; uj{.) — » constant, where the 

constant uj is the neutron detector efficiency. 
Armed with this techique we invert eq. (2): 

i=0 

■w{kn+l)dkn+l ■ ■ ■ w{kn+i)dkn+i (8) 



Where: 



W = —] w — w — \ 



. Taking a simple case, if a; = constant then from eq. (2) 



where n is the measured multiplicity , and with 

1 f 1 

w = 1 = -1 

we can calculate from (8) 

Substituting n + i ^ n we obtain the so called Diven formula Q for 



n>n 



If uj is small we can arrive at a solution containing large oscillating 
and sometimes even negative components of p„ Q] . On the other hand 
it was demonstrated that the moments { < n > and ) of the 
same pn prove to be very stable in the case of multiplicities of fission 
neutrons 0. We show these results in Table 1. 
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Table 1: True moments of distributions for different experiments |Q] 



No. exp. 


M 


No. events 


<n> 




1 


48.3 


7169 


2.690 ±0.036 


1.388 ±0.076 


2 


48.2 


65015 


2.690 ±0.015 


1.290 ±0.025 


3 


44.4 


6928 


2.690 ±0.038 


1.212 ±0.084 


4 


39.9 


20359 


2.690 ±0.025 


1.173 ±0.057 


5 


23.7 


4039 


2.690 ±0.071 


1.230 ±0.272 


6 


22.0 


4039 


2.690 ±0.075 


1.587 ±0.311 



The same conclusion can be drawn from analitical calculation for 
Poisson distribution Q. 

In the general case lj = uj{k) we can prove Q that the true binomial 
moment 

Bj ~ — J wi ■ ■ ■ Wjfj{wi, . . . , Wj)dwi ■ ■ ■ dwj (11) 
where fj is the measured inclusive distribution. 

2 Gamma moments from the data sum- 
mary tape 

u) is the probability of e~^e^ pair creation of a secondary 7: 

u; = 1 - exp (^^^ = ^ (12) 

where is the potential length, L — L{k) is the radiation length and 
w denotes the conversion weight. The general prescription for the true 
binomial moments Q 

Bk - (Bk) (13) 

\ / DST 

where B is the following for every event: 

= (14) 

where a = ^ ^ ^ and n is the detected number of gammas in an 
event and the summation goes for all the a different set of indices. 



4 



E.g. n = 4, 

= 2, a = 6 

B2 = 'Wi'W2 + WiWs + W1W4 + W2W3 + + W3W4 

In addition to Bk we calculated the errors and the correlations of Bk 
from the DST: 

{ABkf = ((Bk-Bk)) (15) 

\ ^ ^ / DST 

ABkABi ^ ({Bk ~ BkKBi - Bi)) (16) 

\ / DST 

Thus we obtained significant resuhs for the first three binomial mo- 
ments. Assuming that all gammas come from neutral pions: 

P2n — Pn I P2JI+1 ~ " 

we can calculate arbitrary 7r° moments using the proper generating 
function: 

G(^)(z) = ^p(7)z" 

and 

G(-")(z) = 5]p(:")z" = G(^)(z^) 

In this way we calculated the average multiplicity in accordance with 

earlier published data § and ^ 2.18+2.81 ^ 2.5 which 

is equal to < n'-'^°'> >= 2.49 ± 0.04 for Tr~p interactions at 40GeV. The 
behaviour of ir^n data on C2 = 1.69 ± 0.11 and C3 = 3.61 ± 0.43 are 
similar to C2 = 1.64 ± 0.07 and C3 = 3.38 ± 0.32 found for n^p data. 

We can compare our results with the 5m hydrogen bubble chamber 
data on ir^p at 250 GeV The statistics is larger (20000 events) but 
< uj >~ 14% is smaller, the experimental detailes have been described 
in§. 

At 250 GeV C2 = 1.55 ± 0.12 and cg = 3.04 ± 0.51. 
Within the errors the KNO moments C2 and C3 do not show a 
violation of KNO scafing ||] between 40 and 250 GeV. 

3 FRITIOF and shifted KW distribution 
for 

We have generated 14500 FRITIOF events at every sample. The 
FRITIOF reproduces the mean multiplicities, the second scaled mo- 
ments C2 and C3, which can be seen in Table 2. 

We should note that a shifted KW distribution has successfuly de- 
scribed ||l^ the KNO moments and the distributions for the single 
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Table 2: Calculated parameters of shifted KW and FRITIOF distributions 

for C3 



Experiment 


1 + e 


A 


< n >* 


C3 [pred] 


C3 [exp] 


C3[FRITI] 


TT-n 40 GeV 


0.955 


0.71 


3.16 


3.57 


3.61 ±0.43 


2.91 ±0.11 


TT-p 40 GeV 


0.965 


0.75 


3.37 


3.38 


3.38 ±0.32 


2.85 ±0.10 


TT-p 250 GeV 


0.978 


0.77 


4.43 


2.98 


3.04 ±0.51 


2.53 ±0.09 



hemisphare data of DELPHI and OPAL collaborations. In order to 
predict the third scaled moments C3-S we use a shifted KW distri- 
bution with parameter m=:2, which has proved to be successful for 
charged particles in inelestic pp collisions ||ll|, |l^ . We carry out a shift 
with +1, since the KW distribution is equal to zero in n=0, and we use 
the so called stick approximation. It means that we use the continuous 
(denoted by *) KW distribution 

Pn = ^^T7TT^^^'"""'exp[-f^z™] 

<n>^r{A) 

where 

r"'(A+i^) n 

F = r-rr^ z ~ m = 2 

r™(A) < n 

Taking the sum of P* for n=0,l,2,... we can define the remaining e in 
the Euler-MacLaurin formula: 



We can form the discrete pobabilities 



Pn = 

fulfilling the requirements: 
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We display the calculated parameters e, < n and A in Table 2. 

Using these parameters we can predict C3, which arc in good agree- 
ment with the true expermental data. With the shifted P*^i we create 
the continuous KNO function 

ij{z) =< n P*+i - (< n > +1)(1 + e)2p„(5, A) 
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where 

n + 1 

< n >*= (< n > +1 1 + e), = z 

A KNO function with parameters A=0.743 and m=2 represents the 
calculated points In Fig.l. KNO scaling is seen as a function of z at 
two energies. 

Conclusions 

1. The use of generating functional technique provides an elegant 
and concise derivation of formulae relating the true distribution 
function to the measured ones. With the aid of this general 
method we have immediately obtained Nifeneckers results on 
the generator functions and Diven |^ formula for the true dis- 
tribution of fission neutron multiplicity, as a special case: w = 
constant. 

2. Since the mean detection efficiency for 7 is 25% in propane (at 
40 GeV) the reconstruction of the true multiplicity distribution 
from the measured one is not efficient, but still the first three 
binomial moments of the true distribution can be obtained. 

3. We have compared our results at 40 GeV with the published 
results at 250 GeV. It was found that the KNO moments C2 and 
C3 of tt" are consistent with KNO scaling within their large (15 
percent) errors up to third moment, in a model-independent way. 

4. The experimental KNO moments are in agreement with the FRITIOF 
simulation and with a shifted KW distribution containing the 
first two binomial moments as an input. It is remarkable that 
this KW distribution can predict the measured C3 . All the points 

of 77° multiplicity distributions calculated as a shifted KW show 
a clear scaling curve between 40 and 250 GeV. 

Acknowledgment 

This work was supported by the Hungarian Academy of Sciences (Grant 
No. OTKA T 016377). 

References 

[1] Kaur and J. D. Swain, in Proceedings of ICHEP, Warsaw - 1996, 
p. 569. 



7 



[2] O. Balea et al, Nucl. Phys. B 63, 114 (1973). 
O. Balea et al, Nucl. Phys. B 83, 363 (1974). 

[3] Brown, Phys. Rev. D 5, 748 (1972). 

[4] L. Diosi, Nucl. lustrum. Methods 138, 241 (1976). 
L. Diosi, Nucl. lustrum. Methods 140, 533 (1977). 

[5] H. Nifoncckcr, Nucl. Instrum. Methods 81, 45 (1970). 

[6] B. C. Diven, Phys. Rev. 101, 1012 (1956). 

[7] R. N. Dakowski et al, Nucl. Instrum. Methods 113, 195 (1973). 

[8] R. N. Diamond et al, Phys. Rev. D 29, 368 (1984). 

[9] Z. Koba, H.B. Nielsen and R Olesen, NPB 40, 317 (1972). 
[10] S. Krasznovszky and I. Wagner, Phys. Lett. B 295, 320 (1992). 
[11] S. Krasznovszky and I. Wagner, Nuovo Cimento A 76, 539 (1983). 
[12] S. Krasznovszky and I. Wagner, Phys. Lett. B 306, 403 (1993). 



8 



